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Abstract This study compares a number of coupled heat 
and mass transfer models and presents numerical compari¬ 
sons of phenomenological coefficients between the four 
models (Stanish, Perre, Pang, and Avramidis) that are most 
frequently used in the literature to describe wood-drying 
processes. The USDA sorption isotherm, the Hailwood- 
Horrobin model, was adopted to calculate the relations 
between moisture content in wood and water vapor pressure 
at any temperature. Due to different assumptions about the 
driving forces of heat and mass transfer, coefficients in each 
model represent different values for moisture content and 
temperature and are closely related to each other. In the 
case of isothermal mass transfer, the moisture diffusion 
coefficient in the transverse directions from the Stanish and 
Pang models increased with decreasing moisture content. 
This contradicts the Avramidis and Perre models and numer¬ 
ous experimental results. Thermal diffusion effects on the 
drying process may not be predominant because the noniso¬ 
thermal state is relatively short. Therefore, the Perre model, 
which does not consider the thermal diffusion effect, has 
been used successfully in the drying simulation. However, it 
may be erroneous in certain cases when the nonisothermal 
state prevails over the system, such as building physics. The 
Pang model cannot explain the phenomena of thermal dif¬ 
fusion and moisture thermodiffusion. It might be reasonable 
to modify the thermal diffusion of the Avramidis model, 
which is lower than that of the Stanish model. The apparent 
heat diffusivity was higher than the true heat diffusivity. 
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Introduction 

Heat and mass transfer modeling is important for improving 
our knowledge of many engineering problems, such as the 
moisture-buffering capacity of wood in building physics, 1 
wood drying, and wood-based panel manufacturing during 
hot pressing. Most of the sophisticated models have been 
developed within the wood-drying research field. Their 
importance is increasing because of the demand for wood 
components as building materials is expected to increase 
in the future and given that well-designed drying process 
control is required at present. 

Moisture in wood exists in three phases: water vapor, 
bound water, and free water above the fiber saturation 
point (FSP). There are two phases of water vapor in the cell 
cavities and bound water in the cell walls, below FSP in the 
hygroscopic range. The maximum bound water content is 
at the FSP and is limited by the number of sorption sites 
available. This may differ among species because their 
chemical compositions vary, but it is usually assumed that 
FSP is 30% water, based on the wood’s dry weight at room 
temperature. 

It is well documented that the driving potential for the 
movement of free water and water vapor is the gradient of 
capillary pressure and vapor pressure, respectively. Fur¬ 
thermore, the flux of water vapor can be divided into the 
convection component of the gas mixture flux and the dif¬ 
fusion component of water vapor under the effect of a gra¬ 
dient in vapor concentration. However, the driving potential 
for the bound water differs according to researchers, and 
can be expressed in terms of a number of different system 
variables. Siau stated that under isothermal conditions it 
does not matter which potential is used in the diffusion 
equation as long as the appropriate diffusion coefficient is 
used. The most commonly used potentials are based on 
the chemical potential and the gradient of bound water 
concentration. 

Under nonisothermal conditions, thermally induced 
mass transfer (thermal diffusion or Soret effect) should be 
taken into account in models. This is because there might 
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be a significant temperature difference between the outer 
and inner surfaces of a building envelope, although this may 
be neglected in wood drying because the temperature 
change duration is relatively small. 

In the literature, we can find three distinct formulations 
of nonisothermal transport models that are widely accepted 
in wood drying and other research fields. First, Avramidis 
et al. 3 developed Luikov-type equations mainly for model¬ 
ing wood drying at low temperatures, and are based on a 
series of studies by Stamm, 4 Choong, 5 and Siau. 2,6 Avramidis 
assumed that mass transfer is driven by moisture and tem¬ 
perature gradients. He used the combined bound water 
diffusion coefficient that involves both bound water and 
water vapor movement. However, the model can be applied 
only below the hygroscopic range and at low temperature 
below the boiling point of water because he did not consider 
capillary movement above the hygroscopic range and water 
vapor convection after change of water state to vapor. 
Second, Perre’s model 7 is based on Whitaker’s multiphase 
transfer model 8 with an additional potential of bound water. 
Perre et al. 7 focused on heat and mass transfer modeling 
above the boiling point of water, in which mass transfer by 
gas convection is dominant. The bound water diffusion 
coefficient is based on Stamm’s result, 4 which is similar to 
the Avramidis model. Third, Stanish 9 assumed that bound 
water diffusion is driven by the chemical potential, which 
makes his model different from Perre’s model. It should be 
noted that both Stanish and Perre developed general models 
for nonhygroscopic as well as hygroscopic materials. Pang 10 
focused mainly on high-temperature wood drying and he 
used the driving potential for bound water similar to the 
chemical potential of Stanish’s model. 

All the above models were compared for their ability to 
predict experimental results. Satisfactory agreements were 
obtained over given material and environmental conditions. 
However, not all the physical and transport coefficients 
were determined experimentally, and the validity would be 
uncontestable. The application of a model to various wood 
species might be questionable under a wide range of dynamic 
external conditions. 

Kamke and Vanek 11 compared the capacities of 12 wood¬ 
drying models with the experimental drying results of mois¬ 
ture change rate and distribution that included the above 
three models (Avramidis, Perre, and Stanish). They stated 
that “none of the models violated any of the known wood¬ 
drying phenomena, but there was a high variability between 


models’ results that was not in close agreement with the test 
runs. It might be due to uncertain coefficients for the models, 
different degrees of simplifications, and different ways of 
solving the heat and mass transfer equations.” In addition, 
the equations for boundary conditions and external heat 
and mass convection coefficients are usually different among 
models, which might affect the numerical results. Yeo and 
Smith 12 developed the convective mass transfer coefficient 
conversion method, which proves that boundary-layer 
theory is useful for evaluating external moisture resistance 
during wood drying, using surface moisture data. 

However, little work has been done on the numerical 
comparison of governing equations for heat and mass trans¬ 
fer models in wood because their assumptions and equation 
forms are different, and it is not easy to compare those dif¬ 
ferences explicitly between models. The objective of this 
study was to numerically compare four models in the hygro¬ 
scopic range after modifying them to the same form of 
partial differential equations. 


Formulation of the heat and mass transfer equation 


Conservation of mass and energy 


Stanish et al. lj and Perre et al. 7 developed the three mass 
conservation equations for the multiphase mass transfer of 
bound water, water vapor, and air, as shown in Table 1 (the 
notation used in the equations and throughout the text 
is defined in Table 2). The fluxes of bound water differ 
between these two reports. In the hygroscopic range, Stanish 
et al. 13 assumed that the flux of bound water is proportional 
to the chemical potential and volume fraction of the cell 
wall, but Perre et al. 7 considered only the bound water 
diffusion. The vapor fluxes can be divided by gas con¬ 
vection and vapor diffusion, and they adopted similar 
assumptions. 

However, the changes in air density and internal gaseous 
pressure can be neglected when modeling heat and mass 
transfer below the boiling point of water. This assumption 
enables the three conservation equations to be simplified to 
two conservation equations and the flux of water vapor can 
be redefined by assuming the air flux is negligible: 13 
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Table 2. Definitions of notation 


Notation 

Definition 

Unit or known value 

Cp 

Specific heat of wood 

J/kgK 

C pa 

Specific heat of air 

J/kgK 

C PV 

Specific heat of water vapor 

J/kgK 

»v 

Bulk binary diffusivity of water vapor in air 

m 2 /s 

D„ 

Diffusion coefficient of wood 

m 2 /s 


Stanish diffusion coefficient 

m 2 /s 

I^bP 

Perre diffusion coefficient 

m 2 /s 


Pang diffusion coefficient 

m 2 /s 

I^bL 

Diffusion coefficient of bound water in longitudinal direction 

m 2 /s 

I^bT 

Diffusion coefficient of bound water in transverse direction 

m 2 /s 

Fa 

Energy to dissociate a mole of sorbed water from sorption sites 

J/mol 

E h 

Activation energy (38500 - 29000 m ) 

J/mol 

E v 

Molar latent heat of evaporation 

J/mol 

E w 

Activation energy of liquid water 

J/mol 

E s 

Molar differential heat of sorption 

J/mol 

^OW 

Specific gravity of oven-dried wood 

Dimensionless 

K 

Enthalpy of bound water 

J/kg 

h. 

Enthalpy of water vapor 

J/kg 

Ah v 

Latent heat of evaporation 

J/kg 

Ah, 

Differential heat of sorption 

J/kg 


Flux of bound water 

kg/m 2 s 

J, 

Flux of water vapor 

kg/m 2 s 

h 

Flux of air 

kg/m 2 s 

K, 

Permeability of water vapor 

m 2 

m 

Fractional moisture content 

Dimensionless 

M w 

Molecular weight of water vapor 

kg/mol 

A 

Partial vapor pressure 

Pa 

p 

vs 

Saturated vapor pressure 

Pa 

R 

Universal gas constant 

8.314 J/mol K 

S V 

Entropy of water vapor 

J/mol K 

T 

Temperature 

K 

a 

Tortuosity, attenuation factor for water vapor diffusivity 

Dimensionless 

V 

Relative humidity 

Dimensionless 

Aff 

Effective thermal conductivity 

W/mK 

0 

Porosity or void fraction 

Dimensionless 

P 

Wood density 

kg/m 3 

Pow 

Wood density of oven-dried wood 

kg/m 3 

Pv 

Density of water vapor 

kg/m 3 

Pw 

Liquid water density 

kg/m 3 

Pa 

Air density 

kg/m 3 

Pvs 

Chemical potential of saturated vapor 

J/mol 

Pb 

Chemical potential of bound water 

J/mol 

Pw 

Chemical potential of liquid water 

J/mol 

Pg 

Dynamic viscosity of air 

Ns/m 2 

Pv 

Dynamic viscosity of water vapor 

Ns/m 2 


Furthermore, the gradient of water vapor can be reduced 
to one of water vapor pressure by using the ideal gas 
law: 7 



The effective vapor diffusivity and bulk binary diffusivity 
of water vapor in air is given by: 


D eff = acj)D v 
D v = 2.20 x 10 


( T Y' 7 Y 101325 Y 


v 273.152 


V 


g 7 


(3) 

(4) 


Bories 14 stated that the attenuation factor for effective 
vapor diffusivity through porous media depends on the 
material’s texture (tortuosity, porosity distribution, anis¬ 
otropy), moisture content, and experimental conditions. 
Akanni and Evans 15 reviewed the tortuosity factor employed 


frequently for porous solids, and represented it as a func¬ 
tion of porosity. However, the coefficient remains unknown 
and differs among researchers. 

Stanish 9 adapted values of a of 0.05 and 0.01 for southern 
pine and Douglas fir, respectively. This study assumed a to 
be 0.03 for the Stanish model. In the Perre model, was 
adapted to be 0.001 to include the effect of porosity, and a 
is much lower than that used by Stanish. 9 

Because this study focuses on comparisons of the models 
only in the hygroscopic range and below the boiling point 
of water, the effect of gas convection can be neglected. 
Therefore, the change rate of water is equal to the total flux 
of bound water and water vapor. The mass transfer equa¬ 
tion can be represented by: 

p B ^ = -Vih + h) (5) 
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The heat-transfer equation can be represented using the 
energy conservation law. 


dT 


pC p — = V ■ (V VT) - V. (J b h b + J A) 


dt 


( 6 ) 


where the density, specific heat, and thermal conductivity 
of wood with changing moisture are calculated using equa- 

-1 S 

tions proposed by Forest Products Laboratory. 
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ow 


Modification of models 


Water vapor pressure is a function of moisture content and 
temperature, and is the product of relative humidity and 
saturated vapor pressure. 


dP dP 

V7 v = —^Vm + — 
dm dT 

P v = (p(m,T)P vs (T) 


( 10 ) 

( 11 ) 


The saturated vapor pressure is a function of tempera¬ 
ture only, but the relative humidity in a material is a func¬ 
tion of both moisture content and temperature, which can 
be obtained only from the measurement of the sorption 
isotherm. Furthermore, vapor pressure with respect to 
temperature is a function of relative humidity and satura¬ 
tion vapor pressure. 


VP v = P vs ^Wm + P vs ^VT + (p^VT 


dm 

VP, = uVm + vVT 


vs 


dT 


dT 


( 12 ) 

(13) 


By substituting Eq. 13 for the gradient of partial 
vapor pressure into the expressions for the fluxes, simulta¬ 
neous heat and mass transfer equations with two state 
variables, moisture content and temperature, can be 
represented. 


'ow 


dm 

dt 

dT 


- V*(flnVm) + V*(% 2 V7) 


pC p — = Via 21 Vm) + V-(a 22 VT) 


(14) 

(15) 


To compare the transfer behaviors between models and 
materials, Eqs. 14 and 15 can be arranged to: 


dm 

~dt 

dT 

dt 


= V-(K u Vm) + V- (K 12 y T) 


= V-(K 21 Vm) + V-(/f 22 Vr) 


(16) 

(17) 


It should be noted that the shrinkage or swelling of wood 
and the change in density and specific heat with time were 
neglected in Eqs. 14 and 15. Equation 17 is not appropriate 


for solving because the density and specific heat are func¬ 
tions of moisture content. 

By some mathematical manipulations, the heat and mass 
transfer coefficients of the Stanish and Perre models can be 
rearranged as shown in Table 3. The four transfer coeffi¬ 
cients given by Pang are available directly from his report, 
and those of Avramidis et al. are simply derived after 
modifying the mass transfer rate of the energy conservation 
equation. 


Results and discussion 

Moisture diffusion 


For the transfer coefficients of the Stanish, Perre, and Pang 
models, the sorption isotherm at a given temperature should 
be known previously. The Hailwood-Horrobin sorption 
model by Simpson 17 was used to calculate the thermody¬ 
namic coefficients ( u , v, 7], £) as shown in Fig. 1. It should 
be noted that the sorption data are approximate means 
of the adsorption and desorption data, which are based on 
spruce grown in North America. Therefore, the model 
may have some deviations for other species, especially in 
hardwoods. 

Changes in 7] and g with temperature are very small, and 
7 ] is always negative over the whole hygroscopic range 
because 




Bound water diffusion within cell wall 

After sealing off the lumens with molten bismuth, Stamm 18 
first measured 7> BL , bound water diffusion coefficients of 
the cell wall in the fiber direction. The adsorption tests 
were conducted at 26.7°C. Using Stamm’s data, Siau 2 con¬ 
verted to D hT by assuming D hT = D hL / 2.5, and obtained the 
Arrhenius-type equation with the least-squares method. 



Fig. 1. Thermodynamic coefficients based on Hailwood-Horrobin 
model by Simpson 17 




Table 3. Heat and mass transfer coefficients below fiber saturation point for each model 
Model K n (m 2 /s) K u (m 2 /s K) 
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The above diffusion coefficients would give lower values 
when applied to wood drying because they are based on 
adsorption data under conditions of low air velocity. 6 In 
addition, the surface resistance was neglected and it was 
assumed that the surface comes to immediate moisture 
equilibrium with the surrounding atmosphere. Further¬ 
more, the molten metal did not fill all the cavities com¬ 
pletely. Nevertheless, this equation has been used widely in 
the field of wood research. 

Perre and Turner 19 also used Stamm’s data to calculate 
the diffusion coefficient of bound water but they did not 
adopt the activation energy, resulting in a slightly different 
form of equation. 


( 


D bT = exp -12.8183933 +10.8951601m- 


v 


4300 T 
T j 



However, the coefficients given by Eqs. 19-21 are close 
to the experimental results found by Stamm 18 at 26.7°C, 
and change little at higher temperatures, as shown in 
Fig. 2. They increase with increasing moisture content and 
temperature. 

The activation energy is one of the kinetic energies that 
is necessary for the bound water to move from one sorption 
site to another, which is lower than the energy needed to 
evaporate bound water. However, its physical meaning 
seems to be unambiguous because it can be measured only 
by diffusion coefficient experiments and their gradient by 
temperature. 



3(lnD) 

3(1//) J 


= RT 


3(lnD)' 
dT 


m 



Nelson 2021 investigated the activation energy of bound 
water theoretically. He stated that it consists of the energy 
needed to dissociate a mole of sorbed water from sorption 
sites ( E a ) and the activation energy of liquid water (£ w ). E a 



Fig. 2. Bound water diffusion coefficient of cell wall in the transverse 
direction with moisture content and temperature 
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consists of E v and RT/n. The former is the energy change 
from the sorbed state to the liquid state, and the latter is 
the energy change from the liquid state to the activated 
state. 

RT RE 

E b = E v + - + E W = E v + —— + 22175 (23) 

n 0.7 

According to Nelson, 19,20 the transverse diffusion coeffi¬ 
cient of bound water can be expressed as: 


7.88 xHT 5 

("0.685 + 0.9m > 

exp 

f-E b ) 

v 2 - 0.9ra/(0.685 + 0.9m) y 

v 2.653 + 0.9/t^ 7 

l RT) 


(24) 


At 26.7°C, Fig. 2 shows that there is good agreement 
between Nelson’s model and Stamm’s experimental results, 
but there is a tendency to underestimate the diffusion coef¬ 
ficient of bound water at lower moisture content levels and 
to overestimate at higher ones. At higher temperatures, 
there is rather good agreement between the three models 
in the range of m > 0.1, but Eq. 24 represents a lower value 
than the other two models when m < 0.1. 

Equation 24 is somewhat different from Eqs. 19-21 
because E v changes exponentially as a function of moisture 
content. The concept of Nelson’s model is useful for esti¬ 
mating E b by wood species because E w may depend on 
the physicochemical properties of wood and a given 
specimen. 


Combined bound water diffusion 


The conductivity of bound water depends on bound water 
diffusion through the cell walls and water vapor diffusion 
into the cell cavities. As discussed by Choong, 5 the trans¬ 
verse moisture movement mechanism is dominated by 
cross-wall diffusion, which increases as moisture contents 
increase; whereas the longitudinal mechanism is dominated 
by vapor diffusion in the cell lumens, which decreases as 
moisture contents increase. This is supported by many 
experimental results, as shown in Fig. 3 and Table 4, and 
also by Danvind and Ekevad. 22 

The theoretical diffusion coefficient from the Avramidis 
model used in this study originates from Siau. 2 It can be 
predicted more precisely for softwood, such as Scotch pine 
than for hardwood, such as aspen and red oak. 

For convenience, the specific gravity of wood is assumed 
to be 0.40 for the calculation of heat and mass transfer coef¬ 
ficients in this study. The combined or effective transverse 
diffusion coefficients of bound water (K n ) for each model 
are shown in Fig. 4. The parameters necessary for the 
models are adopted in the literature as: 


D bS (kg-s/m 3 ) = 3.0xl0“ 13 

f 4300 'N 

D bP (m 2 /s) - expl -9.9 + 9.8m- 

D bR (kg-s/m 3 ) = 8.0 x 10“ 13 and K v (m 2 ) 


Stanish model 
Perre model 
= 1.0 x 10” 15 


Pang model 


The transverse effective moisture diffusion coefficients 
in the Avramidis and Perre models are similar and increase 



Moisture content 


Fig. 3. Comparison of transverse effective diffusion coefficients by 
theory and experiments during desorption at 60°C. Filled circles, Scots 
pinep filed triangles, Scots pine; open triangles, Scots pine; filled 
squares, aspen (transverse); 28 open squares, red oak (transverse) f filled 
diamonds, theoretical (SG 0.30); 2 open diamonds, theoretical (SG 
0.65) 2 



Fig. 4. Transverse effective moisture diffusion coefficients of wood 
with oven drying and specific gravity of 0.4 at 20°C and 80°C 


with increasing moisture content, because it is assumed that 
Perre included the porosity effect on £> bP , and that the water 
vapor effect through the pit, D v , is very small because = 
0.001. K n in the Avramidis model is based on Siau, which 
did not consider the effect of tortuosity and assumed that 
a - 1 . 0 . 

However, the coefficients from Stanish’s and Pang’s 
models decrease with increasing moisture content because 
D hS and D hR are constant, but £, u , and D w decrease with 
increasing moisture content, which disagrees with most of 
the previous experimental results. All the models consid¬ 
ered in this study assume that all the phases are at the same 
temperature and that the partial vapor pressure is equal to 
its equilibrium pressure, whose condition is said to be in 
local thermodynamic equilibrium. 

Stamm 4 proposed an electrical model to simulate mois¬ 
ture transfer in the hygroscopic range (Fig. 5), shown in 
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R1 [] R3 j R4 



Fig. 5. a Primary paths for tangential and radial moisture transport 
through a cell wall and b secondary paths through a pit 5 




(d) 92 9i 

—Wv— wv 



Fig. 6a-g. The electrical analog models for the paths of tangential 
moisture transport 


Fig. 6a. Figure 6b shows a simplified model integrating sec¬ 
ondary paths through a pit. Because the resistance of T3 

r\ 

and T4 is large, Siau neglected their effects and adopted 
the model shown in Fig. 6d. As Siau stated, the effect of T3 
may be important in woods with high permeability and high 
specific gravity, particularly at low moisture content and 
high temperature. The contribution of the pits in hardwood 
is expected to be less than in softwood because of the 
smaller pit openings in hardwoods. 

However, Perre et al. 7 considered T3 in their model and 
used the model from Fig. 6e, although Fig. 6f may be more 
suitable for woods according to Stamm’s theory. Stanish 
and Pang considered the movement of bound water and 
water vapor in the cavity by including the contribution of 
the pits as parallel paths. They adopted the model in Fig. 
6g, and water vapor transfer through both the cavity and pit 
is more dominant than bound water. Therefore, their models 
are subject to some limitations for wood in the transverse 
direction and may be more suitable for moisture transfer 
through wood in the fiber direction, wood-based panels 
(e.g., medium-density fiberboard and particleboard), and 
homogeneous porous material. For the models of Stanish 
and Pang, the dependence of combined diffusion coefficient 
on moisture content does not change even if the bound 
water and water vapor transfer is transformed to series dif¬ 
fusion such as Fig. 6d. 



Fig. 7. Thermal diffusion coefficients of wood with oven-dried specific 
gravity of 0.4 at 20°C and 80°C 


Thermal diffusion coefficient 

Moisture is transported in the direction from higher tem¬ 
perature to lower temperature, although there is no mois¬ 
ture gradient within a material, and is known as thermal 
diffusion or the Soret effect. Thermal diffusion is more 
important in building physics than drying because wood in 
service is rarely maintained in an isothermal state, and the 
nonisothermal state is relatively short during drying. 

It should be noted that the Avramidis model originates 

r\ 

from Siau, and Perre’s model did not consider the thermal 
diffusion effect. Stanish et al. compared their model to the 
experimental data from the literature. 

The thermal diffusion coefficients ( K n ) of the three 
models are shown in Fig. 7, except for Pang’s model, whose 
value is several orders of magnitude lower than K n . These 
coefficients are found to increase with increasing moisture 
content and are observed to be higher in higher tempera¬ 
ture ranges. Pang’s model has a negative value in the entire 
hygroscopic range because -D bR { 1 - (j))r\ > vK v p v lju v , which 
contradicts the phenomenological observation that mois¬ 
ture flows in the direction of decreasing temperature. The 
Stanish model also has a negative value at low moisture 
content due to the large entropy value of water vapor. 

In the steady state, the Soret effect can be calculated by 
using dm/dT - -K u /K n (Fig. 8). In the nonisothermal state, 
the moisture transfer rate decreases with increasing -dm/dT 
during drying, when the environmental temperature is 
higher than the material temperature. 

The Avramidis model shows a higher value than Stan- 
ish’s model at 20°C, but it is close to the Stanish model at 
80°C because there is little change with temperature in the 
Avramidis model. These two models show that the value 
increases with moisture content; however, it decreases with 
moisture content in the Perre model, and gives a constant 
value near zero over the whole moisture range. Therefore, 
Perre’s model cannot explain the Soret effect in steady state 
experiments. 












































































274 


Figure 9 shows the results of comparison of the Avrami- 
dis, Stanish, and other models, and the relevant equations 
from the literature are shown in detail in Table 5. 

Stanish 9 compared his model to Models 2 and 7 in the 
unsteady state. As he stated, the accuracy of Model 7 must 
be questionable because it has a negative value over the 
whole moisture range. Peralta and Skaar 2j compared the 
steady-state uniaxial moisture profiles to those predicted by 
five different theoretical models. Models 9 and 5 provided 
the best results. Model 4 gave the next best prediction, fol¬ 
lowed by Model 6 (which is based on the Stanish model), 
and the last was Model 2. Siau and Avramidis 24 compared 
three models (2, 8, and 9); Model 9 provided the best fit to 
the experimental results, and Model 2 gave the poorest 
results. 

There is little difference between Models 5 and 9 below 
moisture content of 0.2%, but they deviate greatly above 
that. Model 1 does not consider the material’s hygroscopic 
properties, so with hygroscopic material -dm/dT should be 
lower than that in Model 1. Avramidis and Hatzikiriakos 25 
compared Models 4 and 9 with drying experiments and 
concluded that the thermodiffusion coefficient of Model 9 
resulted in the best fit with the moisture and temperature 



Fig. 8. Plot of -dm/dT against moisture content with oven-dried spe¬ 
cific gravity of 0.4 at 20°C and 80°C 


data. Therefore, it may be concluded from the above 
research that Model 9 is the best one to use for the noniso- 
thermal modeling. 

Equation 23 and Model 5 show that the diffusion coeffi¬ 
cient of bound water decreases, and the Soret effect 
increases as the differential heat of sorption increases. 

Heat diffusion and moisture thermodiffusion 

The same values of volume heat capacity ( pC p ) and heat 
transfer coefficient (A e// ) were adapted to all the models. As 
shown in Table 3, however, the additional term of heat dif¬ 
fusion coefficient ( K 22 ) made the thermal conductivity vary 
between models because the condensation or sorption heat 
occurs in the direction of increasing temperature as the 
moisture content increases or decreases. As shown in Fig. 
10, Perre’s and Pang’s models at all temperatures, and Stan- 
ish’s model at 20°C, show similar trends, with K 22 decreasing 
at increasing moisture content, which was very close to the 
changes in X eff /pC p . However, in Stanish’s model it remains 
constant with moisture content at higher temperature 
(80°C), whereas in the Avramidis model there is a trend of 


0.08 



I I I I I I 
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Moisture content 

Fig. 9. Plot of -dm/dT against moisture content with oven-dried spe¬ 
cific gravity of 0.4 at 50°C 


Table 4. Comparison of parameters in effective diffusion equation D b = D 0 exp(aX + b) from experimental results in the literature 


D 0 (m 2 /s) 

a 

b 

Experimental conditions 


Literature cited 




Direction 

Temperature (°C) 

Species 


325 x 1CT 10 

-28.2 

0 

Longitudinal 

30 

Western fir 

Choong 5 

171 x 10'“ 

-13.3 

0 

Longitudinal 

23.9 

Yellow poplar 

Rosen 30 

0.53 x lO -10 

0 

0 

Tangential 

20 

American beech 

Skaar 31 

0.27 x 10'“ 

6.5 

0 

Radial 

20 

American beech 

Skaar 31 

0.044 x 10'“ 

13.8 

0 

Tangential 

26.7 

Sitka spruce 

Stamm 4 

0.27 x 10'“ 

18.4 

0 

Tangential 

40 

Yellow poplar 

Comstock 12 

5.55 x 10' 4 

16.8 

-5280 IT 

Transverse 

43 

Aspen 

Simpson and Liu 28 

19.2 x 10' 4 

1.15 

-5280 IT 

Transverse 

43 

Red oak 

Simpson 29 

16.0 x 10' 4 

1.45 

-5280 IT 

Transverse 3 

43 

Red oak 

Simpson 29 

6.34 x 10' 14 

5.46 

2.54 x 10” 2 T 

Transverse 

20/50/80 

Scots pine 

Hukka 26 

4.06 x 10- 14 

5.32 

2.66 x 10” 2 T 

Transverse 

20/50/80 

Norway spruce 

Hukka 26 


a Adsorption 
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Table 5. Comparison of the ratio of moisture content gradient to temperature gradient ( dm/dT) from the literature 


Model 

Equation 


Literature cited 

1 

V 



Bramhall 33 

2 

u 

E„ 


Skaar and Siau’ 4 


RT 2 IdEA 



3 

i 

V 

n \ dm Jr 

( T ^ + E s + RTln(<p) ] 

\ dT ) 


Siau 35 

4 

V 

C 

( r dT +£ ‘ + £v + £s + Rnn(,p) ) 

Siau and Zin l6 

5 


c 

E v + E s -C pv T 


Nelson 20 ’ 21 

6 

(-1.038 - 0.049E) (-RT\ncp)T 
<pTS v 

Peralta and Skaar 23 

7 

i 

V 

C 

^T^ + Eb + Es + RT Info)) 

Siau 37 

8 

V 

c 

[ T ^ + E b + RT\n{<p) 

\ dT / 


Siau and Avramidis 24 

9 

c 

<pE b 


Siau and Avramidis, 24 Peralta and Skaar 23 


c 




C — RT 2 


fdqr\ 

\dmj 


T 



Fig. 10. Heat diffusion coefficient of wood with oven-dried specific 
gravity of 0.4 at 20°C and 80°C. Filled circles , Avramidis model; open 
circles , Perre model \ filled triangles , Stanish model; open triangles , Pang 
model; crosses , Eq. 25 


increasing K 22 at higher moisture contents. If K 22 from Eq. 
25 were adopted, the heat diffusion coefficient would exhibit 
a lower value than in the Avramidis model 



1 L + EbpCp_Eb 

pC p { eff 0.018C p T 


KOP 

RT 


^ dm V 

U<pJ, 



The Dufour effect, or moisture thermodiffusion, which 
causes the heat flux resulting from moisture diffusion, is 
shown in Fig. 11. However, there is little experimental evi¬ 
dence of the Dufour effect because it cannot be easily mea¬ 


Fig. 11. Moisture thermodiffusion coefficient of wood with oven-dried 
specific gravity of 0.4 at 20°C and 80°C 


sured as the heat diffuses so much more rapidly through 
wood than does moisture. Therefore, it was difficult to 
determine what model was more appropriate for explaining 
moisture thermodiffusion. In addition, the effect of mois¬ 
ture thermodiffusion on the heat transfer rate may not be 
significant because the heat diffusion is much faster than the 
moisture diffusion. 

The moisture thermodiffusion coefficient ( K 21 ) from the 
Avramidis model increased with increasing moisture 
content, which is due to a function of the activation energy 
of bound water. For the other three models, it decreased 
with increasing moisture content because the enthalpy of 
bound water ( h b ) increases with moisture content. 
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Conclusions 

The driving potential for bound water may differ among 
researchers, and can be expressed in terms of a number of 
different state variables. Under isothermal conditions, it 
does not matter which potential is used in the diffusion 
equation as long as the appropriate diffusion coefficient is 
used. 

Except for the Avramidis model, the other three models 
(Stanish, Perre, and Pang) were developed mainly for mod¬ 
eling high-temperature drying above the FSP The transfor¬ 
mation of the four models with different state variables to 
the same forms results in different phenomenological trans¬ 
port coefficients by moisture content and temperature. 
There is a certain degree of confusion about the material 
parameters. This confusion arises from the use of different 
mathematical expressions and the driving forces assumed. 
In the case of the isothermal mass transfer coefficient, the 
moisture diffusion coefficient in the transverse direction by 
the Stanish and Pang models increased with deceasing 
moisture content. This contradicts the Avramidis and Perre 
models as well as numerous experimental results. 

The effect of thermal diffusion on the drying process 
may not be predominant because the nonisothermal state 
is relatively short. Therefore, Perre’s model, which does not 
consider the thermal diffusion effect, has been used success¬ 
fully in the drying simulation. However, it may be errone¬ 
ous in certain cases when the nonisothermal state prevails 
over the system, such as building physics. The thermodiffu¬ 
sion coefficient, such as from Model 9, may be involved in 
Perre’s model. 

Model 9 is more suitable than Model 4 for the Avramidis 
model. The water vapor movement through the pit should 
be involved in case of wood with high permeability and high 
specific gravity, particularly at low moisture content and 
high temperature, when calculating the moisture diffusion 
coefficient (Fig. 6b). 

Therefore, the application of a model to various wood 
species would be questionable under a wide range of exter¬ 
nal conditions, although all four models were verified by 
experimental results, and satisfactory agreements were 
obtained for a given material and set of environmental 
conditions. When developing and adopting new models, 
examination of whether the transport coefficients have 
physical meaning and are appropriate over the range of 
moisture content and temperature under consideration is 
required. 
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